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Abstract 

We present a global fit of current available experimental results on Xc decays 
within next-to-leading-order perturbative QCD. The quality and reduced errors of 
recent data improve the agreement between theory and experiment. 
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The study of charmonium has recently received renewed interest. All aspects, from 
spectroscopy, to decays and production have benefited from new data and novel theoret- 
ical developments. Recent experimental data concern measurements of total |T], ^[ and 
radiative decay widths ||, §, |5|, the observation of the 1 P\ state ||, measurements of 
decay branching ratios of B mesons into P-wave charmonium states J7J , production cross 
sections in hadronic collisions, both at fixed target [§] and at collider energies || pi], [T| . 
New theoretical developments cover, among other things, attempts to extract the value 
of a s from lattice QCD calculations of quarkonium decay widths [0, the development 
of a systematic approach to the problem of infrared ambiguities in the production and 
decay of P-wave states |13|, [1^] , discovery of new production mechanisms at high p t JT5 . 



The large disagreement observed between expected and measured production rates at 
the Tevatron [JTU|, has also led to new speculations about the existence of exotic spec- 



troscopy in the charmonium system fl6f . We believe that this discrepancy provides an 



important arena in which to test our understanding of the boundary domain between 
perturbative and non-perturbative QCD. In this respect, it is fundamental to derive the 
best possible perturbative predictions, in order to be able to firmly assess the need for, 
and to properly model, possible additional non-perturbative phenomena required for a 
complete description of this physics. With this goal in mind, we recently completed a 
full next-to-leading-order (NLO) calculation for the total production cross section of 3 Po,2 
states in hadronic collisions. A consistent use of this calculation, requires the inclusion 
of phenomenological parameters extracted from a NLO analysis of inclusive decay widths 
of these states. In this letter, we therefore present an updated comparison between QCD 
predictions for the decay widths of 3 Pj (x c j) states and the latest experimental data. We 
will not include the 1 Pi state here, as the available data are still insufficient. 



Several detailed studies of this subject have appeared in the past [17], [13], [18], [L9| . We 
feel that the new data justify an update. Furthermore, we improve the commonly used 
expression for the decay width of the xa state to light hadrons |20H , with the inclusion of 
a finite 0(a 3 ) term which has usually been neglected. 

We will show that the inclusion of all available data, in addition to relaxing some theory 
constraints used in the literature so far, allows a consistent global fit in terms of three 
parameters. The first parameter is simply a s , the QCD coupling constant. The second 
parameter, |i?'(0)| 2 , corresponds to the derivative of the non-relativistic P-wave function 
at the origin. The third parameter is required to regulate the infrared divergency which 
appears in the standard calculation of P-wave decays to a gqq final state. This parameter 
is an infrared cutoff, usually loosely referred to as the binding energy E^i^. It appears in 
the standard expression for the decay widths f2(| as a factor L = log[4m^/(4m 2 — M 2 )] = 
log(M/2Ei,i n d), where m c is the constituent charm quark mass. In the recent formulation 
of P-wave decays introduced in ref. JT3I , this logarithm is absorbed into the color octet, 3 Si 



component of the cc wave function, which turns out to be proportional to the combination 
a s \R'(0)\ 2 L. As discussed in [14], this relation provides a solid basis for the study of higher 



order perturbative corrections, as well as providing a rigorous framework to regulate IR 
divergencies appearing in the evaluation of the qq^gxj cross section. 



We shall start by collecting here the expressions for NLO x decay widths [20] that will 
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be used in our fit: 

r(xj-ry) = A7$aZj^W(l + BA) (J = 0,2) (1) 

fil r \ 7T / 



r( Xj ^LiJ) = A^ffi^fl + Cj^l+Tr^fli, (J = 0,2) (2) 



. , 56rif 9 a s |-R'(0)| 2 



c 



where: 



A? = 27 AV = f Bb = y-f ^ = "f (4) 
Ag 9 = 6 Af = | C = 8.772 C 2 = -4.827 n f = 3 (5) 



To follow the spirit Ref. |T3|, we introduced the parameter ifg; which we define as: 



if. = (6) 

The precise connection between this parameter and the color octet wave function, can be 
found in [[L3|]. 



Following the suggestion of JUj, we chose to scale the widths by the charm constituent 
quark mass, m c , rather than by the mass of the quarkonium states. This is consistent 
with the neglect of higher order non-relativistic (NR) corrections, and with the inclusion 
of all non-perturbative effects into |i?'(0)| 2 . Since |_R'(0)| 2 and m c always appear in the 
fixed combination \R'(0)\ 2 /m^, the precise value of the charm mass will not change the 
fitted value of the other two free parameters, a s and H%. It is important to point out, 
furthermore, that it is possible to make consistent use of the extracted value of |_R'(0)| 2 
in the calculation of production cross sections, by properly identifying the origin of the 
mass terms present in the theoretical cross sections. For example, LO direct P-wave 
production cross sections are proportional to \R'(0)\ 2 / M 7 . Of the seven powers of mass, 
four have the same origin as those appearing in the decay widths. The remaining three 
are from phase space. It is therefore consistent to write the overall factor appearing in 
production as |-R'(0)| 2 /16m^/M 3 , where M is the charmonium state mass. With this 
choice, no ambiguity in the choice of m c arises when using our fitted values in the study 
of production cross sections. 

In our expressions we neglected the very small contribution of the 3 gluon decay of 
the XcitHI- We included however a finite, non-logarithmic contribution to T(xci^LH), 
which was not evaluated in the original calculation ||20| . To our knowledge this piece has 



always been neglected in previous studies |. It turns out that this term is not numerically 

2 We thank E. Braaten for pointing out this fact and suggesting we evaluate and include this 
contribution. 



2 



negligible when compared to the formally dominant logarithm, and should therefore be 
included. 

Corrections to the NR approximation are not known, but are expected to at least 
partly cancel when taking ratios of widths. In a recent series of papers [|L8|| , Consoli and 
Field argued that by properly taking into account the phase space reduction due to the 
effective mass of the gluon, it is possible to provide a consistent perturbative description 
of the decay widths of charmonium and bottomonium S-wave states, without need for the 
inclusion of significant NR corrections. This is consistent with the observation that EM 
decays (such as r\ c — ^77 or ip— >£ + £~) are correctly predicted if the 0{v 2 ) corrections are 
set to 0. It is not known whether this applies to P-waves as well, but we will comment 
on the consequences of these ideas for our fits of x decays at the end. 

Taking ratios or differences of appropriate widths, can lead to expressions which are 
independent of |-R'(0)| 2 or H 8 , or both. These ratios can be used to estimate with smaller 
theoretical uncertainty the only really perturbative parameter of the theory, namely a s . 
Some examples which are often used in the literature are: 

r(Xco~>77) /15\ 1 + m 



T( Xc0 ^LH)-T(xci^LH) /15\ 1 + (8.772 + 28/27)^ 



r(xc 2 - 


■*LH) - 


r(xd- 






r(xco- 


^77) 




r(xc2- 




r(x c i- 






r(x C 2- 


^77) 




r(xc2- 




r(xci- 





4 / 1 + (-4.827 + 35/9)^ 



/135A ,a 

V 8 ) QK a s ' 1 + (-4.827 + 35/9)^ 1 ; 



2 J Q a s ' 1 + (-4.827 + 35/9)^ v ; 

Notice that the inclusion of the non-logarithmic contribution to T{xc\— >LH) (the factor 
35/9 appearing in the above equations) significantly reduces the C(af) corrections to the 
difference r(x c2 — >LH)-T(xci^LH). 

In reference [|13| the value of a s was extracted from the measurement of bottomonium 
decays, evolved down to the charm mass. This implicitly assumes that a s has to be 
evaluated at the scale of the heavy quark mass, and led to a value of a s (m c )=0.25±0.02 
at m c =1.5 GeV. We prefer instead to leave a s as a free parameter, to be fit together with 
|-R'(0)| 2 /m^ and Hg. In fact we believe that the safest and less restrictive assumption on 
the scale of a s is that it has to be same for all decay processes, therefore enabling us to 
use the same value of a s regardless of J and of the final state. This choice is also free of 
ambiguities related to the actual value of m c . While this will not allow us to extract a 
value of Aqcd, h provides however a less restrictive constraint on the comparison of data 
with QCD. We will verify at the end that the fitted value of a s is in reasonable agreement 
with what expected fixing the renormalization scale to be of the order of the charm mass. 
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Another choice performed in JT3| was to use only the leading order (LO) version of 
equations (f7|)— (]ToD . This was justified by the fact that the contributions to decay widths 
via the color-octet 3 Si component of the wave function are only known to LO. However, 
higher order corrections to the terms proportional to H$ do not depend on J, as they arise 
from the S'-wave component of the wave function [T4\. We can absorb these universal 
higher order corrections into the parameter H$, so that their effect does not change the 
results of the fit for the variables a s and |i?'(0)| 2 . We therefore feel that it is justified to 
use the full 0(a z s ) expressions given above. 

We collect the data that will be used in our fit in Table 1. We notice that a key 
measurements, namely r(x c0 — ^77), was never actually published pif . The large error 
associated to it will not weigh this measurement significantly in our fit, which is unfortu- 
nate since the ratio of the photonic widths of the even x states is a very good probe of the 
theory. The other debated item is r(x c2 — ^77), whose central value differs significantly 
among various experiments || [|, 0] (we do not quote previous older results, some of which 
are simply upper limits). It is not our duty to judge on the value of the experimental 
data, so we chose, contrary to other authors JIB] , to include them all and let the associated 
errors drive the fit. 

We provide here the results of the fit to these seven measurements, using the theoretical 
widths provided in equations (ffl)-(H): 



a s = 0.286 ±0.031 ^ ^ = 0.60 ± 0.10 MeV H 8 = 4.2 ± 0.7 MeV (11) 

16m* 



E x = 13.7 ± 2.3 MeV xju = 7A ( 12 ) 

The total x 2 °f the fit is 7.1 for the four degrees of freedom. The poor quality of the fit 
is mostly due to the discrepancy between the different measurements of r(x c2 — ^77). We 
included, for reference, also the numerical value of H±, defined in [[13] as: 

* = I ™ (13) 

We notice that the value of ot s returned by the fit is consistent with a s {m c ) = 0.30, 
obtained at two loops using m c =1.5 GeV and the value of A4~ /oop = 235 MeV extracted 
from DIS data [p2fl . The value of the derivative of the wave function is also consistent 
with potential model calculations (see for example the recent update in ref . |£| ) . 

In Table 2 we provide the distribution of differences between expected and measured 
decay widths, relative to the experimental errors, for the seven measurements considered. 
The theoretical values are obtained using the fitted values of parameters. It is interesting 
to see what happens if one of the two best measurements of T(x C 2— is removed from 
the fit. Removing the CLEO data point gives as central values for the fit : 

\R'(0) I 2 

a, = 0.298 ± 0.034 ' 1 '} = 0.55 ±0.11 MeV H 8 = 3.9 ± 0.7 MeV (14) 

16m* 
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#1 = 12.6 ± 2.3 MeV x% = 3-5 



(15) 



The central values have not changed significantly from the global fit, but the x 2 is now 
consistent with the 3 remaining degrees of freedom. The CLEO measurement in this case 
would be off by less than 2 sigma from the theoretical expectation. 

Removing the E760 data point gives as central values for the fit : 

\R'(0) I 2 

a s = 0.195 ±0.031 - — ^ = 1.23 ± 0.3 MeV H 8 = 7.7 ±2.1 MeV (16) 

16m^ 



Hi = 28.3 ± 7.8 MeV xju = 31 ( 17 ) 

The central values have now moved significantly. The x 2 has improved, thanks to the 
larger relative error quoted by CLEO for r(x c2 — ^77). Notice that in this case the value 
of a s extracted is consistent with what determined by CLEO in their analysis of their 
measurement (a s =0.219±0.127 ||). The E760 measurement is this case would be off 
by about 6 sigma from the theoretical expectation based on the values of parameters 
extracted from the fit. 

In Table 3 we present the comparison between the width ratios given in equations (|7])- 
(|10|) and the data. The experimental error bars are quite large, due to the propagation 
of errors in ratios of differences. It is likely that some of the systematics or statistical 
errors are correlated and will cancel in these combinations, but we did not pursue this 
possibility in absence of enough details on the experimental analyses. The same results, 
derived by excluding either E760 or CLEO from the fit, are also included in Table 3. 

As a final exercise, we include in our analysis the effect of an effective gluon mass 
on the final state phase space. Following ref. |IB| , we applied a correction factor to the 
hadronic decay widths. Notice that contrary to 3 S\ decays, where the final state involves 
three gluons at LO, there are only two gluons at LO in xo,2 decays. As a consequence, the 
impact of this correction is less significant. It is not clear to us what is the right procedure 
to extend this idea to decays to qqg final states. Since these are dominated by the soft 
gluon region, where the quarks carry most of the energy, and since this domain is already 
screened by the IR cutoff, we chose not to include any correction factor for these final 
states. We collect here the results of the global fit, which are only meant to be indicative 
of the possible size of these effects: 

\B'(Q) I 2 

a s = 0.326 ±0.030 J — ^ = 0.69 ± 0.11 MeV H 8 = 4.2 ± 0.5 MeV (18) 

16m^ 



H x = 15.8 ± 2.0 MeV xju = 6 - 6 ( 19 ) 

The most significant change is in the value of a s , as already pointed out in ref. fT8|| . 
Notice also a slight improvement in the quality of the fit. We also mention that a clear 
prediction of the Consoli and Field approach is that the ratio of hadronic widths of xo and 



5 



X2 should be insensitive to the effective gluon mass. The experimental value of this ratio 
is 7.9 ± 3.3. The results of our fit from Eq. (|i~T|) yield a theoretical ratio of 5.4 ± 1.3, those 
from Eq.(|TE|) yield 6.4 ± 1.4. Both results are in agreement with the data, the second one 
being slightly better because of the larger value of a s . 

In conclusion, we find that current data on 3 Pj charmonium decays are well consistent 
with NLO perturbative QCD. The error bars are still large for more incisive tests of the 
theory, and leave room for deviations from the naive NR approximation and for higher 
order perturbative corrections. Nevertheless the agreement found is encouraging, and 
hopefully relieves serious concerns raised in earlier works. A reduction in the significant 
discrepancy found between the r(x c2 — ^77) widths measured by different experiments will 
be of fundamental importance to guide the extraction of theoretical parameters from the 
data. Likewise, a new measurement of r(x c0 — ^77) would be very helpful. 

Acknowledgements: We are very grateful to E. Braaten and M. Consoli for several 
discussions and precious comments to this work. 
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Process 


T (MeV) 


Reference 




T(Xco- 




(4.0 ± 2.8) x KT 3 


Crystal Ball 


21 




r(xc 2 - 


m) 


(0.321 ±0.095) x 10-3 


E760 § 




r(xc 2 - 


m) 


(1.08 ±0.38) x 10- 3 


CLEO § 




r(xc2- 


m) 


(3.4 ± 1.9) x 10-3 


TPC2 7 [§ 




r(xco- 


>LH) 


13.5 ±5.4 


Crystal Ball 


21 


1 


nxd- 


>LH) 


0.64 ±0.10 


E760 Q 




r(xc 2 - 


>LH) 


1.71 ±0.21 


E760 





Table 1: Most recent experimental results on Xc decay widths. The errors were obtained 
by combining in quadrature the statistical and systematic errors given in the quoted 
references. The widths to light hadrons were obtained from total widths by removing the 
contributions of known radiative decays. 



Process 


Dat a-Theory /Error 


r(x c0 ^77) 


-0.44 


T(x c2 ^77) (E760) 


0.55 


r(x c2 ^77) (CLEO) 


-1.8 


T(x c2 ^77) (TPC2) 


-1.6 


T(xco^LH) 


-0.74 


t( Xc1 ^lh) 


-0.13 


T(xc2^LH) 


0.29 



Table 2: Fractional differences, relative to the experimental errors, between data and 
theory predictions after the global fit. 
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Theory 


Data 






Global 


No E760 No CLEO 




r(Xc0^77) 
r(Xc2^77) 




7.4 ±0.8 


5.7 ±0.45 7.7 ±0.9 


12.5 ± 9.5 (E760) 
3.7 ±2.9 (CLEO) 
1.2 ± 1.0 (TPC2) 


r( Xc0 ^L/j)-r( Xcl - 


>LH) 


7.7 ±0.5 


6.4 ±0.4 7.9 ±0.5 


12.0 ±5.7 


T{ Xc2 ^LHyT{ Xcl - 


>LH) 


r(Xco^77)xio 3 




2.4 ±0.5 


5.0 ±1.5 2.2 ±0.5 


3.7 ±2.7 


T{ Xc2 ^LH)-T{ Xcl - 


>LH) 


r(Xc2^77)x!o 4 




3.3 ± 1.0 


8.8 ±3.4 2.9 ±1.0 


3.0 ±1.1 (E760) 
10 ±4 (CLEO) 
32 ± 19 (TPC2) 


r( Xc2 ^L/7)-r( Xcl - 


>LH) 



Table 3: Comparison between data and theory for ratios of widths (includes global fit, fit 
without E760 and without CLEO r(x c2 — ^77) datum). 
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